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Arbuscular mycorrhiza (AM) is the most wide-spread plant-fungus symbiosis on earth. Investigating 
this kind of symbiosis is considered one of the most promising ways to develop methods to nurture 
plants in more natural manners, avoiding the complex chemical productions used nowadays to pro- 
duce artificial fertilizers. In previous work we used the Calculus of Wrapped Compartments (CWC) 
to investigate different phases of the AM symbiosis. In this paper, we continue this line of research 
by modelling the colonisation of the plant root cells by the fungal hyphae spreading in the soil. This 
study requires the description of some spatial interaction. Although CWC has no explicit feature 
modelling a spatial geometry, the compartment labelling feature can be effectively exploited to de- 
fine a discrete surface topology outlining the relevant sectors which determine the spatial properties 
of the system under consideration. Different situations and interesting spatial properties can be mod- 
elled and analysed in such a lightweight framework (which has not an explicit notion of geometry 
with coordinates and spatial metrics), thus exploiting the existing CWC simulation tool. 

1 Introduction 

Arbuscular mycorrhiza (AM) is the most wide-spread plant-fungus symbiosis on earth [25 ]. Investigating 
this kind of symbiosis is considered one of the most promising ways to develop methods to nurture plants 
in more natural manners, avoiding the complex chemical productions used nowadays to produce artificial 
fertilizers. 

In previous works lfT71l45l in the context of the BioBITS project iflTTl we investigated different phases 
of the AM symbiosis by using the Calculus of Wrapped Compartments (CWC) |fl9l , a variant of the 
Calculus of Looping Sequences (CLS) |7l|6|. These calculi have been designed to represent in a formal 
way biological entities such as molecules, DNA strands and cells. Starting from an alphabet of atomic 
elements {atoms for short) and from an alphabet of labels, CWC terms are defined as multisets of atoms 
and compartments. A compartment consists of a wrap (a multiset of atoms detailing the elements of 
interest on the membrane), a content (a term) and a type (a label). The evolution of the system is driven 
by a set of rewrite rules modelling the reactions of interest (elements interaction and movement). In 
previous work, we have introduced in CWC rules a notion of rate which allow CWC to be equipped 
with a semantics given in terms of Continuous Time Markov Chains (CTMCs), thus enabling stochastic 
simulation built following the standard Gillespie's approach 11241 . Moreover, we have defined a hybrid 
simulation algorithm to speed up the simulation procedure while keeping the exactness of Gillespie's 
method, |[T8l . The simulation tool [2] has been integrated in the FastFlow @ programming framework 
for the execution of parallel simulations in multi-core platforms Q . 

*This research is funded by the BioBITs Project (Converging Technologies 2007, area: Biotechnology-ICT), Regione 
Piemonte. 
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In this paper, we continue our ongoing research on the AM symbiosis by modelling the colonization 
of the plant root cells by the fungal hyphae spreading in the soil. The interaction begins with a molecular 
dialogue between the plant and the fungus. Host roots release signalling molecules characterized as 
strigolactones which induce alterations in the fungal physiology and hyphal branching. In response to 
the plant signal AM fungi produce a symbiotic signal, the "Myc Factor". One of the first response 
observed in epidermal plant root cells is a repeated oscillation of Ca 2+ concentration preparing to the 
reorganization of the cell for the fungal penetration. Following this event, fungal hyphae enter and 
cross the epidermal cells, growing inter-and intracellularly all along the root in order to spread fungal 
structures. Once inside the inner layers of the cortical cells, specialised fungal hyphae differentiate into 
arbuscules where the nutrients exchange between the two organisms occur. 

To represent this kind of behaviour we need to model some spatial properties of the system like the 
movement of the hyphae in the soil and their penetration and articulation in the root tissue. To this aim 
we exploit the fact that the notion of compartment can also be used, in a very natural way, to represent 
spatial regions (with a fixed topology) in which the labels play a key role in defining the spatial properties. 
These kind of compartments are thus considered as the topological sectors in the analysed model. In this 
framework, the movement and growth of system elements are described via specific rules (involving 
adjacent compartments) and the functionalities of biological components can be affected by the spatial 
constraints given by the sector in which they interact with other elements. A similar approach can be 
found in f34l where the topological structure of the components is expressed via explicit links which 
require ad-hoc rules to represent movements of biological entities. 

Although our approach does not support an explicit representation of continuous space, we can ex- 
ploit the topological structure of compartments, their inter-relations, and their relative placement as a 
starting point to implicitly define spatial properties. This level of abstraction is suitable to model some 
situations where it is not needed to explicitly specify the space occupied by the elements of the sys- 
tem. Avoiding to introduce explicit features for modeling a spatial geometry (e.g., coordinates, position, 
extension, motion direction and speed, rotation, collision and overlap detection, communication range, 
etc.) may in some cases reduce the complexity of the models and simplify the analysis of the system and 
the simulation procedure. 

1.1 Summary 

The remainder of the paper is organised as follows. In Section [2] we recall the syntax and semantics 
of CWC. In Section [3] through a few chosen paradigmatic examples, we start giving some hint about 
spatial modelling and analysis within the CWC framework. In Section [4] we exploit the approach to 
model some spatial interaction in the Arbuscular Mycorrhizal Symbiosis. In particular, we show how 
molecules (carrying the communication signals between the plant and the fungus) diffuse in the soil, and 
how the fungus contacts, through its hyphae, the cells of the plant root. Finally, in Section [5] we briefly 
discuss related and possible future work and draw our conclusions. 

2 The Calculus of Wrapped Compartments 

The Calculus of Wrapped Compartments (CWC) (see |[T9l[T7l ) is based on a nested structure of ambi- 
ents delimited by membranes with specific proprieties. Biological entities like cells, bacteria and their 
interactions can be easily described in CWC. 
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Figure 1: (a) represents (a b c\ (b) represents (a b c\(d e\»Y'Y; (c) represents 
(abc\(de\*Y fg) ( 

2.1 Syntax 

Let be a set of atomic elements {atoms for short), ranged over by a, b, and J2? a set of compartment 

types represented as labels ranged over by £, A term of CWC is a multiset t of simple terms where a 

simple term is either an atom a or a compartment ia\ t'Y consisting of a wrap (represented by a multiset 
of atoms aj^ represented by a (possibly empty) and a content, represented by a term t') and a type, 
represented by the label I. An empty multiset is represented symbolically with 

As usual, the notation n * t denotes n occurrences of the simple term t. An example of term is 
2* a b (c d\e fp representing a multiset (multisets are denoted by listing the elements separated by 
a space) consisting of two occurrences of a, one occurrence of b (e.g. three molecules) and an £-type 
compartment (c d\e f) e which, in turn, consists of a wrap (a membrane) with two atoms c and d (e.g. 
two proteins) on its surface, and containing the atoms e (e.g. a molecule) and / (e.g. a DNA strand). See 
Figure [T] for some other examples with a simple graphical representation. 



2.2 Rewriting Rules 

System transformations are defined by rewriting rules written over an extended set of terms including 
variables. The l.h.s. component of a rewrite rule p is called pattern and the r.h.s. component o of a 
rewrite rule is called open term. Patterns and open terms are multiset of simple patterns p and simple 
open terms o defined by the following syntax: 



P 
o 



a 
a 



(a x\ p X) 
(a x I oY 



X 



where a is a multiset of atoms, x and X are variables that can be instantiated, respectively by a 
multiset of atoms (on a wrap) or by a term (in the compartment contents). The label I is called the 
type of the pattern. Patterns are intended to match, via substitution of variables, with simple terms and 
compartments occurring as subterms of the term representing the whole system. Note that we force 
exactly one variable to occur in each compartment content and wrap. This prevents ambiguities in the 
instantiations needed to match a given compartment. |^] 

A rewrite rule is a triple (£,p,o, denoted by £ : p 1 — > o), where p and o are such that the variables 
occurring in o are a subset of the variables occurring in p. The application of a rule £ : p 1 — > o to a term 



'the wrap of a compartment contains the elements characterizing the membrane enclosing it. So we assume that it cannot 
contain other compartments but only atomic elements (for instance proteins). 

2 The linearity condition, in biological terms, corresponds to the exclusion of transformations depending on the presence of 
two (or more) identical (and generic) components in different compartments (see also [39|). 
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t is performed in the following way: 1) Find in t (if it exists) a compartment of type I with content u 
and a substitution a of variables by ground terms (i.e. terms without occurrences of variables) such that 
u = o(p X), where X is a variable not occurring in p, o; and 2) Replace in t the subterm u with o(o X). 
We write 1 1 — > t' if t' is obtained by applying a rewrite rule to t. 

For instance, the rewrite rule £ : a b i — > c means that in all compartments of type £ containing an 
occurrence of a b, it can be replaced by c. While the rewrite rule £ : (x\a b X) i[ i — > (x\c X) i2 means 
that, if contained in a compartment of type £, all compartments of type £\ and containing an a b in their 
content can change their type to £2 and an occurrence of a b can be replaced by c. 
Remark. For uniformity we assume that the term representing the whole system is always a single 
compartment labelled T with an empty wrap, i.e., all systems are represented by a term of the shape 
(• J f) , which we will also write as t for simplicity. 

2.3 Stochastic Simulation 

A stochastic simulation model for biological systems can be defined along the lines of the one presented 
by Gillespie in l24ll . which is, de facto, the standard way to model quantitative aspects of biological 
systems. The basic idea of Gillespie's algorithm is that a rate function is associated with each considered 
chemical reaction which is used as the parameter of an exponential distribution modelling the probability 
that the reaction takes place. In the standard approach this reaction rate is obtained by multiplying the 
kinetic constant of the reaction by the number of possible combinations of reactants that may occur in 
the region in which the reaction takes place, thus modelling the law of mass action. In |fl9ll , the reaction 
rate is defined in a more general way by associating to each reduction rule a function which can also 
define rates based on different principles as, for instance, the Michaelis-Menten nonlinear kinetics. 

In the examples of this paper it is reasonable to assume a linear dependence of the reaction rate on 
the number of possible combinations of interacting components, so we will follow the standard approach 
in defining reaction rates. Each reduction rule is then enriched by the kinetic constant k of the reaction 
that it represents (notation £ : p 1 — > o). For instance in evaluating the application rate of the stochastic 

rewrite rule R = i : a b X 1 — > c X to the term t = a a b b in a compartment of type £ we must consider 
the number of the possible combinations of reactants of the form a bint. Since each occurrence of a can 
react with each occurrence of b, this number is 4. So the application rate of R is k ■ 4. 

2.4 The CWC simulator 

The CWC simulator is a tool under development at the Computer science Department of Turin Uni- 
versity, based on Gillespie's direct method algorithm [24]. It treats CWC models with different rating 
semantics (law of mass action, Michaelis-Menten kinetics, Hill equation) and it can run independent 
stochastic simulations over CWC models, featuring deep parallel optimizations for multi-core platforms 
on the top of FastFlow [4]. It also performs online analysis by a modular statistical framework. 

3 Topological Interpretations within CWC 

In this section we provide some hint about spatial modelling and analysis within the CWC framework 
by means of some paradigmatic examples. For simplicity, we will consider the compartments in our 
system to be well-stirred. Thus, also the compartments representing the spatial sectors of our topology 
will consist of well-mixed components. Note that the topological analysis described in the following 
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examples could be expressed also with other calculi with compartmentalisation such as, for example, 
BioAmbients [42], Brane Calculi (H) and Beta-Binders EQ. 

3.1 Cell Growth and Proliferation: A CWC Grid 

Eukaryotic cells reproduce themselves by cell cycle which can be divided in two brief periods: the inter- 
phase during which the cell grows, accumulating nutrients and duplicating its DNA and the mitosis (M) 
phase, during which the cell splits itself into two distinct cells, often called "daughter cells". Interphase 
proceeds in three stages, Gl, S, and G2. Gl (Gapl) is also called the growth phase, and is marked by 
synthesis of various enzymes that are required in 5 phase, mainly those needed for DNA replication. In 
S phase (Synthesis) DNA replication occurs. The cell then enters the G2 (Gap2) phase, which lasts until 
the cell enters mitosis. Again, significant biosynthesis occurs during this phase, mainly involving the 
production of microtubules, which are required during the process of mitosis. 

Spatial representation of the growth of cell populations could be useful in many situations (in the 
study of bacteria colonies, in developmental biology, and in the analysis of cancer cells proliferation). 

Intuitively, we may express cell proliferation by describing a finite space through a CWC grid such 
that on each grid cell we can insert a cell of the growing population under analysis. An k x n grid could 
be modelled in CWC with k * n compartment labels modelling the grid cells. A cell may grow in one of 
the adjacent grid cells (with a given scheme of adjacency modelled by the set of rewrite rules). 

Intuitively, using the atomic element e to denote an empty grid cell, an empty grid could be modelled 
by the following CWC term: 

Prol_Grid= (Je) 1 ' 1 ... (Je) 1 '" 
(j,) 2 ' 1 ... {\ef- n 

{\e) mA ... (\e) mM 

If we represent a cell of our population as a compartment with membrane m and with a single atom 
in its content representing the phase of the cell cycle (M, Gl, S, G2), a grid cell i,j may contain the cell 
(m\M) ce " that could grow, with kinetic k m i t , towards the empty grid cell i, j + 1 with the rewrite rulej^j 

T : (x\ (m y\ M Y) cM X) ! ^(zJ e Z)^ +l (xj (m y\ Gl Y) ce " X)^(z\ (mj Gl) ceU Z)^ +1 
while cells may change their internal state (cell cycle phase) with the rules: 

ij : (m x\ Gl X) ceU (m x\ S X) cdl 
i,j : (m x\ S XY M fe 2 (m x\ G2 X) cM 
ij : (mx\G2X) ce " k ^ (mx\M X) ceH 

Note that each element of the grid can be occupied only by one cell, so this model allows only a 
limited growth of the cells, corresponding to the dimension of the grid. However this could be enough to 
investigate the properties of the cell growth. 



note that under our assumptions the variables will always match only with the empty multiset 
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Figure 2: Quorum Sensing in a distributed topology. 



3.2 Quorum Sensing and Molecular Diffusion 

Molecular diffusion (or simply diffusion) arises form the motion of particles as a function of temperature, 
viscosity of the medium and the size (or mass) of the particles. Usually, diffusion explains the flux of 
molecules from a region of higher concentration to one of lower concentration, but it may also occur 
when there is no concentration gradient. The result of the diffusion process is a gradual mixing of the 
involved particles. In the case of uniform diffusion (given a uniform temperature and absent external 
forces acting on the particles) the process will eventually result in a complete mixing. Equilibrium is 
reached when the concentrations of the diffusing molecules between two compartments becomes equal. 

Diffusion is of fundamental importance in many disciplines of physics, chemistry, and biology. In 
this subsection we will consider the diffusion of the auto-inducer communication signals in a quorum 
sensing process where the bacteria initiating the process are distributed in different sectors of an environ- 
ment (see Figure [2]). We consider the sectors of the environment to be connected through channels with 
different capacities and different speeds regulating the movement of the auto-inducer molecules between 
the different sectors. 

The CWC model describing the system can be found in Appendix A. We modelled the four sectors as 
four different compartments with special labels Si allowing the flow of the oxo3 auto-inducer molecules 
with different speed. Bacteria (modelled as compartments with label m) were marked as "active" when 
they reached a sufficient high level of the oxo3 auto-inducer. The results of a simulation running for 30 
seconds is shown in Figure [3] 

Although CWC does not have a direct mechanism to formulate the spatial position of its terms and a 
metric distance between them, it has the potential to describe topologically different spatial concepts. For 
example, the compartments are able, through the containment specification, to express discrete distance 
among the entities that are inside and those on the outside. Moreover, through the rule kinetics we are 
able to model space anisotropy (lack of uniformity in all orientations), which are useful to formulate 
complex concepts of reachability. The model discussed in this section is a clear example of these issues. 
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Figure 3: Simulation results of oxo3 molecules (left figure) and active bacteria percentage with respect 
to the total number of bacteria (right figure) in the different sectors. 

4 Spatial Interactions in the Arbuscular Mycorrhizal Symbiosis 

Arbuscular mycorrhizae, formed between 80% of land plants and arbuscular mycorrhizal (AM) fungi 
belonging to the Glomeromycota phylum, are the most common and widespread symbiosis on our planet 
l25l . AM fungi are obligate symbionts which supply the host with essential nutrients such as phos- 
phate, nitrate and other minerals from the soil. In return, AM fungi receive carbohydrates derived from 
photosynthesis in the host. AM symbiosis also confers resistance to the plant against pathogens and 
environmental stresses. Despite the central importance of AM symbiosis in both agriculture and natural 
ecosystems, the mechanisms for the formation of a functional symbiosis between plants and AM fungi 
are largely unknown. 



4.1 The Biological Interactions under Investigation 

The interaction begins with a molecular dialogue between the plant and the fungus ll25ll . Host roots 
release signalling molecules characterized as strigolactones CD which are sesquiterpenes derived from 
the carotenoid pathway li33~l l32ll Within just a few hours, strigolactones at subnanomolar concentrations 
induce alterations in fungal physiology and mitochondrial activity and extensive hyphal branching ifTOl . 
Root exudates from plants grown under phosphate-limited conditions are more active than those from 
plants with sufficient phosphate nutrition, suggesting that the production of these active exudates in roots 
is regulated by phosphate availability ll36l l47l . Strigolactones also appear to act as chemo-attractants: 
fungal hyphae of Glomus mosseae exhibited chemotropic growth towards roots at a distance of at least 
910 /im in response to host-derived signals, possibly strigolactones ll44l . Some solid evidence has been 
presented for AM fungal production of a long-hypothesized symbiotic signal, the "Myc Factor", in re- 
sponse to the plant symbiotic signal. Fungal hyphae of the genus Gigaspora growing in the vicinity of 
host roots, but separated from the roots by a membrane, release a diffusible substance that induces the 
expression of a symbiosis-specific gene in Medicago truncatula roots ll29l . This expression was corre- 
lated both spatially and temporally with the appearance of hyphal branching, and was not observed when 
hyphal branching was absent. The external signal released by these fungi is perceived by a receptor on 
the plant plasma membrane and is transduced into the cell with the activation of a symbiotic signalling 
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pathway that lead to the colonization process. One of the first response observed in epidermal cells of the 
root is a repeated oscillation of Ca 2+ concentration in the nucleus and perinuclear cytoplasm through the 
alternate activity of Ca 2+ channels and transporters ll38ll30l[T6l . Kosuta and collegues ||30l assessed my- 
corrhizal induced calcium changes in Medicago truncatula plants transformed with the calcium reporter 
cameleon. Calcium oscillations were observed in a number of cells in close proximity to plant/fungal 
contact points. Very recently Chabaud and colleagues [16] using root organ cultures of both Medicago 
truncatula and Daucus carota observed Ca 2+ spiking in AM-responsive zone of the root treated with 
AM spore exudate. 

Once reached the root surface the AM fungus differentiates a hyphopodium via which it enters the 
root. A reorganization of the plant cell occurs before the fungal penetration: a thick cytoplasmic bridge 
(called pre-penetration apparatus -PPA-) is formed and it constitutes a trans-cellular tunnel in which the 
fungal hypha will grow [23 ]. Following this event, a hyphal peg is produced which enters and crosses the 
epidermal cell, avoiding a direct contact between fungal wall and host cytoplasm thanks to a surrounding 
membrane of host origin ll22l . Than the fungus overcomes the epidermal layer and it grows inter-and 
intracellularly all along the root in order to spread fungal structures. Once inside the inner layers of the 
cortical cells the differentiation of specialized, highly branched intracellular hyphae called arbuscules 
occur. Arbuscules are considered the major site for nutrients exchange between the two organisms and 
they form by repeated dichotomous hyphal branching and grow until they fill the cortical cell. They are 
ephemeral structures ll46l l27ll : at some point after maturity, the arbuscules collapse, degenerate, and die 
and the plant cell regains its previous organization |[T2l . The mechanisms underlying arbuscule turnover 
are unknown but do not involve plant cell death. In arbuscule-containing cells the two partners come in 
intimate physical contact and adapt their metabolisms to allow reciprocal benefits. 

4.2 The CWC Model 

We model the Arbuscular Mycorrhyzal symbiosis in a 2D space. The soil environment is partitioned into 
different concentric layers to account for the distance between the plant root cells and the fungal hyphae 
where the inner layer contains the plant root cells. The plant root tissue is also modelled as concentric 
layers each one composed by sectors. 

For the sake of simplicity we modelled the soil into 5 different concentric layers (modelled as differ- 
ent compartments labelled by L0-L4 identifying L4 with the top level T). Although the root has a complex 
layered structure we simplified it into 3 concentric layers (one epidermal layer and 2 cortical layers) each 
one composed by 5 sectors mapped as different compartments labelled by e, (i = 1 ... 5) for epidermal 
cells and Cjt (j = 1,2 and k = 1 ... 5) for cortical cells. Exploiting the flexibility of CWC, compartment 
labels are here used to characterize both their spatial position and their biological properties. All the 
compartments and the main atoms of the model are depicted in Figure |4| 

The stigolactones (atoms S) derived from the carotenoid pathway (atom Cp) are released by epidermal 
cells and can diffuse through the soil degrading their activity when farthest from the plant. This situation 
was modelled by the following rules: 

(Ras) Lq : (x\CpX) e -^S(x\CpXY- i = 1...5 

(R s .) L i+l : (*J S X) L ' A S{x\ X) L - i = 0. . . 3 (1) 

(R S4 ) T : sA. 

Since phosphate (atoms P) conditions regulate strigolactones biosynthesis inhibiting the carotenoid 
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Figure 4: 2D spatial model of the Arbuscular Mycorrhyzal symbiosis. 



pathway (atom ICp), the following rules were considered : 



(A 



ICPtj 



P(x\ Cp X) e ' P(x\ IC P X) e - i = 1 . . .5 



(2) 



The fungal structure is described by atoms where each atom represents a particular type of fungal 
component. The atom Spore represents the fungal spore while atoms Hyp represent fungal hyphae. AM 
fungal hyphae branching towards the plant root is regulated by the strigolacones activity. If we locate 
the fungal spore at a fixed distance from the plant root (let us say at the layer where < k < 4) then 
the spore hyphal branching is generated stochastically in the relative inferior and superior layers. In our 
experiments the fungal spore was placed at the soil layer Li- The hyphal branching is radial with respect 
to the spore position causing a branching away from the spore and promoted at the proximity of the root 
by the strigolactones. This branching process was represented by the following rules: 



(R Hj ) 



L 2 

T 

u 



Spore(x\ X) 



Spore(x\ Hyp X) L < 



(x\ Spore X) L2 ^ Hyp(x\ Spore X) L * 
(xj Hyp X) L i A Hyp(x\ Hyp X) Li 



(3) 



S Hyp(x\ X) Lo 1-4- S Hyp(x\ Hyp X)^> 



The AM fungal production of Myc factor (atom Myc) at the proximity of the root cells inducing the 
calcium spiking phenomenon (atom Co) at a responsive state (atom Resp) which prepares the epidermal 
cells to include the fungal hyphae was modelled by the following rules: 



{RMyc) Lo 

{RMycCcij ) Lo 
{Rcai) U) 



Hyp i — — ^ Myc Hyp 

Myc(Respx\ X) e > (Resp x\ Ca X) e - i=\...5 
Hyp(Resp x\ Ca X) e ' ^ Hyp(NResp x\ Hyp X) e ' i = 1 . . .5 



(4) 
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where the epidermal cells change their state from responsive to not responsive (atom NResp) since 
they can be penetrated by only one hypha. 

Once the fungal hyphae penetrate the epidermal cells they are allowed to branch in the two directions 
of the underneath cortical layer and from this cortical layer to the inner one. The cortical cells change 
their state from responsive (atom Resp) to transitive (atom Trans) and then to mycorrhized (atom Arb) 
since an hypha can grow in the cortical intercellular spaces even if already mycorrhized but only one 
arbuscule can be generated. The hyphal branching from the epidermal layer to the cortical layers is 
modelled through the following rules: 



(R CBu ) Lq : (NResp x\ Hyp X) e '(Resp y\ Y) c " ^> (NResp x\ Hyp X) e '(Trans y\ Hyp Y) c ^ 
i = 1 . . . 5,7 = i(mod5),i+ \(mod5); 

(R C Bi 2l ) U) : (Trans x\ Hyp X) Ci > (Resp y\ Y) C2 > ^> (Trans x\ Hyp X) c " (Trans x\ Hyp Y) c v 
i = 1 . . .5,j = i(mod5),i + \(mod5)\ 

(RcB2 2 ) L : (Arb x\ Hyp X) c "(Resp y\ Y) c v ^ (Arb x\ Hyp X) Cl '(Trans y\ Hyp Y) c v 
i = 1 . . .5,7 = i(mod5),i+ \(mod5); 

(5) 

Finally, the modelling of each fungal hyphae penetrated into a cortical cell to create an arbuscule 
(atom A) and start the symbiosis is performed by the following rules: 

(R AiJ ) L : (Trans x\ Hyp X) c -J ^ (Arb x\ Hyp X) C 'J i = l,2;j = 1...5 (6) 

4.3 Simulation Results 

The initial term describing our system is given by T: 

T= (j (j (\Spore(\ (j T ')^) L ^y 
where 
T' = V U 

V = V l V 2 V 3 V 4 V 5 (7) 
Vi = (Resp\ Cp) e - 

U = U n U l2 U n U u U l5 U 2 i U 22 U 23 U 24 U 25 
U jk = (Resp\ 



The model has been used to virtually simulate 21 days of the fungus-plant interactions. The reaction 
rates have been chosen to resemble the experimental observations of "in vivo" fungal germination, hy- 
phal branching and arbuscules creation under similar environmental conditions of the simulations. The 
rates were set as following: K C p = 0.8, K s = 0.1, K ICP = 0.01, K F = 0.2, K H = 0.05, K Hl 0.0l, K M = 0.2, 
Kmc = 0.2, Kc a = 1-0, Kqb = 0.9, Ka = 1.0. In Figure[5]we show a surface plot representing an exempli- 
fication of the CWC simulation about strigolactones diffusion (on the left) and fungal hyphal branching 
(on the right) along the soil layers during time. 

We compared our simulation results adding the phosphate (100 P atoms) to the soil layers yielding 
to the results given by Figure [6] In this case we have a lower strigolactones production by the root cells 
and an hyphal branching not promoted in the direction of the root. 

We performed 60 simulations in the three conditions of low phosphate (10 P atoms), medium phos- 
phate (50 P atoms) and high phosphate (100 P atoms) conditions leading to a 60% of arbuscules in the 
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Figure 5: Simulation results of strigolactones diffusion (left figure) and fungal hyphal branching (right 
figure) along the soil layers during time. 
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Figure 6: Simulation results of strigolactones diffusion (left figure) and fungal hyphal branching (right 
figure) along the soil layers at high phosphate concentrations. 

first case, 30% in the second case, and 10% in the latter case. The first arbuscules occur around the 
8 th — 9 th day of simulation at low phosphate condition while they are formed from the 9 th to the \5 th day 
at high phosphate conditions. These results reflect the in vivo experiments on several plants inoculated 
with different AM fungi 1551 151 Egl. 

Our long term goal is to design a virtual spatial model that would enable us to test "in silico" various 
hypotheses about the interaction between the physiological processes that drive the arbuscular mycorrizal 
symbiosis and the signalling processes that take place in this symbiosis. 

5 Conclusions, Related and Further Work 

In the last few years, computer scientists have developed several languages and tools to build and anal- 
yse models of biological systems, often relating different behaviour that operates on different levels of 
abstraction and on different time and spatial scales HTl l43l l20l [T4l I2TI . 



14 



Modelling Spatial Interactions in the AM Symbiosis using CWC 



This work constitutes a first step towards providing CWC with the features required by the emerging 
field of spatial systems biology, which aims at integrating the analysis of biological systems with spatial 
properties. 

For the well-mixed chemical systems (even divided into nested compartments) often found in cellular 
biology, interaction and distribution analysis are sufficient to study the system's behaviour. However, 
there are many other situations, like in cell growth and developmental biology where dynamic spatial 
arrangements of cells determines fundamental functionalities, where a spatial analysis becomes essential. 
Thus, a realistic modelling of many cellular phenomena requires space to be taken into account |[28l . 

There are several cases in which the spatial properties of a system could be analysed by just focusing 
on a finite (and quite limited) number of sectors. If we focus ourselves in this kind of situations, CWC 
appears to be expressive enough for quite effective descriptions. Indeed, the considerations presented 
in this paper about the topological organisation of the paradigmatic examples shown in Section [3j also 
hold for other calculi which are able to express compartmentalisation (see, e.g., BioAmbients H2l . Brane 
Calculi 031, Beta-Binders EQ, etc.). 

Von Neumann' Cellular Automata ll37l were developed as a computational tool inspired by biological 
behaviour and later used to model biological systems. Cellular Automata are defined by a finite grid of 
cells, with the state of each cell, evolving in discrete time, affected locally by just the state of the nearby 
cells. We can mimic the framework of Cellular Automata within CWC by defining a grid in a way similar 
to what we have done in Section [3TTI 

There are other situations, however, where a geometrical distance becomes necessary beyond a purely 
topological organization (for example in a more detailed analysis in the field of developmental biology). 

As a consequence, the need of an extended spatial analysis in systems biology has brought to the 
extension with (even continuous) spatial features of many formalisms developed for the analysis of bi- 
ological systems: see, e.g., geometric process algebra 031, spatial P-sy stems 0, spatial CLS and 
Bioshape |fl3l . 

Theoretically, building a more complex structure over the approach presented in this paper, we could 
define a geometric space, with coordinates and distance metrics. We might define, e.g., on the lines of the 
approach presented in [34], a surface language modelling a 2D grid or a 3D space over the CWC model. 
While adding too many features to the model (e.g., coordinates, position, extension, motion direction 
and speed, rotation, collision and overlap detection, communication range, etc.) could rise heavily the 
complexity of this kind of approach, a detailed study on a subset of these features could still be carried 
out on several systems. We plan to investigate this line of development in near future work. 
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Appendix A: Quorum Sensing and Molecular Diffusion Model 

Pseudomonas aeruginosa uses quorum sensing to keep low the expression of virulence factors until the 
colony has reached a certain density, when an autoinduced production of virulence factors is started. The 
bacterium membrane, denoted as m, contains only a DNA strand. The DNA is modelled as a sequence of 
genes lasO lasR lasl (atom LasORI), where lasO represents the target to which a complex autoinducer/R- 
protein binds to promote transcription. The following rewrite rules model the system behaviour (the 
kinetics are taken from [6]): 



(*i) 
(*3) 
(*5) 

(Ri) 

(*n) 

(*13) 
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oxo3 LasR i— > R3 
R3 LasORI R03 
R03 ^ R03 LasR 
Lasl i — y • 
oxo3 i — > • 



(8) 



Rules R\ and R2 describe the production from the DNA of proteins LasR and Lasl, respectively. Rule 
R3 describes the production of the autoinducer, modelled as atom oxo3, performed by the Lasl enzyme. 
Rules R4 and R$ describe the complexation and decomplexation of the autoinducer and the LasR protein, 
where the complex is denoted R3. Rules R-j, R%, Rg describe the binding of the activated autoinducer 
(the R03 complex) to the DNA and its influence in the production of LasR and Lasl. Finally, rules 
Rio, Ru, R\2 describe the degradation of proteins. In particular, rule 7?i2 models the degradation of the 
autoinducer, which can happen both inside and outside the bacterium. Bacteria are marked as "active" 
(atom A) when reaching a sufficient high level of auto-inducer inside the respective sector otherwise they 
are marked as "inactive" (atom I). 

The following rules describe the ability of the autoinducer to cross the membrane, in both directions 
inside the sectors (different kinetics are used to model the different bandwidths of the channels). An 
autoinducer inside the bacterium is modelled as a non-positional element, while autoinducers outside 
have an associated spatial sector S; for i = 1, ... ,4: 
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ND 
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ND 
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ND 



Table 1 : Values for the kinetic rates fey 



The diffusion of the autoinducer across the different sectors is described by the following rules: 

T : (xj oxo3 X) s -(y\ Y) S J ^% (xj X) s '{y\ oxo3 Y) s j (10) 

where kij is related to the sectors channel capacities. In our experiments the values of kij are shown 
in Table [T] where fey = fey,. 

The simulations were carried out using the following initial term T (see Figure [2]): 

7i = (J 10* (J/ LasORI) m ) Sl T 2 = (J 15* (J / LasORI) m ) Sl 

T 3 = (J 15* (J/ LasORI) m ) Si T 4 = (J 8 * ( J / LasORI)"') S4 (11) 

r = 7i r 2 r 3 r 4 

The model above can be easily enriched by adding rules describing also the movement of the bacteria 
between the different sectors and bacteria replication. 



